Milne phase for the Coulomb quantum problem related to 
Riemann's hypothesis 
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' Abstract. We use the Milne phase function in the continuum part of the spectrum of the 

particular Coulomb problem that has been employed by Bhaduri, Khare, and Law as an 
equivalent physical way for calculating the density of zeros of the Riemann's function on 
the critical line. The Milne function seems to be a promising approximate method to calculate 
^ , the density of prime numbers. 
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, From a 1995 PRE paper of Bhaduri, Khare, and Law [1] one can obtain the following formula 
G^ ■ for the density of zeros of Riemann's zeta function on the critical line 
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q: -^(^) = -^ + ^1^4* (7 + ^0 )1' (1) 
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P_i, where the digamma function is the logderivative of the gamma function, ^(-z) = F (z) /T{z). 

Formula 16 of the same paper gives the phase shift of a repulsive Coulomb potential obtained 
from an inverted oscillator with a hard wall at the origin for the unconventional value of the 
partial wave number I = —\- Taking the derivative of that formula, one obtains 
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nc{e) = ^ + T^Re ^ - + 2- , (2) 



F(e) = — — tan~"^(cosech7re). (3) 

Thus, under an appropriate shift, the two expressions given by (1) and (2) differ only by an 
exponentially small term as noted by the Indian team of authors. 

Our aim in this work is to apply another technique based on the so-called Milne phase 
function for the calculation of the 'density of states' in the continuum of the same Coulomb 
problem. Indeed, being a phase, one might think a priori that Milne's function has something 
to do with the nontrivial zeros of the Riemann function. Not only this is a different procedure, 
but it might be a quite competitive approximation for nz- Previously, Korsch [2] applied the 
same method with very good results in the case of bound states for a few illustrative cases 
of elementary quantum mechanics. He also established that in this approach the density of 
quantum states is nowhere unique and not even necessarily positive (!?). For the Coloumb 
problem under focus here the nonuniqueness issue is actually an advantage because in a 
certain sense one can choose by trial and error a better Milne approximation for nz- The 
Milne function is expressed through the following formula 

nM{y, e) = ^ = ^ -2 , (4) 
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where p is the solution of the Pinney nonlinear equation and 0i and 02 are the two linear 
independent solutions of the repulsive Coulomb problem, i.e. 



ky-- \n{2ky) - ^ + ArgF (^Z, + 1 + 

and 



02(?/, e) = cos f A;?/ - I H2ky) - ^ + ArgF {lr + l + |) 
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where, as we mentioned before, = —j. The other symbols are as follows: e = — ^ is a 
reduced spectral parameter, where E is the spectral parameter in the initial inverted oscillator 
problem and to is the angular frequency of the oscillator; k = y — x"^, where x > is 
the oscillator coordinate. The superposition constants a and P are determined by arbitrary 
'initial' conditions at some point t/o- We fix them in the most 'economic' way, which is at the 
point yo = ^ allowing one to eliminate the logarithm in the argument of the trigonometric 
functions. Moreover, we employ the Eliezer-Gray prescription, for details see the master 
thesis of Espinoza [3]. Thus 

-*4)-ta(i + ^A.grQ + |)), (7) 

A three-dimensional plot of nM{y,£) is displayed in Fig. 1, where an expected oscillatory 
behaviour can be seen. Not only the procedure based on the Milne phase function could 
compete very well with other approximate methods for the density of prime numbers but 
there is a further advantage on which we briefly comment in the following. As is well known, 
the Milne phase function enters as a basic ingredient in the Ermakov-Lewis phase- amplitude 
approach for parametric oscillator problems (for recent applications see [4]). To transform the 
Coulomb problem at hand 
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dy"^ L y 16y^ ^ ^ ^ 

into a parametric dynamical problem for a unit mass classical particle, one can use the well- 
known map to canonical classical variables (f) — q and ^ — P leading to 

«^^ = P (10) 

dy 

dy L y Iby^-' 

where the coordinate y plays the role of the classical Hamiltonian time. Various quantities, 
such as the Ermakov-Lewis invariant and geometrical angles can be calculated easily for this 
particularly interesting Coulomb problem (due to its connection with prime numbers). But 
we want to emphasize a different point here. In the parametric oscillator interpretation and of 
y as Hamiltonian time, one can use nuiy,^) as a direct tool for time-energy (and therefore 
time-imaginary axis) characterization of the density of prime numbers as can be inferred from 
Fig. 1. 
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Figure 1. nniy, e) = l/p^ for y G [0.1, 10] and e e [0.1, 10]. 
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